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I.      ABSTRACT 


In  this  study,  the  flow  characteristics  of  air  around  a 
circular  cylinder  were  determined  by  two  means:  by 
analytical  solution  of  the  Na vier-Stokes  equations,  and 
actual  measurement  of  the  flow  itself.  A  DISA  LD  5586  laser 
docplsr  anemometer  system  was  employed  to  measure  normal  (v) 
and  parallel  (u)  velocity  components  of  a  laminar  flow 
regime  about  a  three-sixteenths  inch  diameter  cylinder 
mounted  in  a  32  x  45  inch  wind  tunnel  (Reynolds  number 
approximately  40)  . 

A   finite  element  grid  was  constructed  using  triangular 

elements  which  encompassed  the  upper   half-cylinder.    Flow 

measurements   using   the   laser  velocimeter  were  confined  to 
this  region  cf  interest. 

Theory  was  compared  against  experimental  results,  and 
the  feasibility  of  the  utilization  of  this  particular  type 
of  laser  anemometer  was  evaluated. 
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II.      DESCRIPTION  OF  LASER  ANEMOMETER 

A  basic  laser  anemometer  is  comprised  of  four  primary 
elements: 

(1)  a  coherent  light  source 

(2)  optics 

(3)  signal  converter  (photomultiplier  tube) 

(4)  signal  processor/display 

The  principle  of  the  laser  anemometer  is  to  create  a 
grid  pattern  and  measure  the  speed  with  which  a  particle  in 
a  fluid  free  stream  passes  through  the  grid.  The  following 
assumptions  were  taken  into  consideration: 

(1)  particles  injected  into  the   free  stream  are 

small  enough  to  faithfully  follow  stream 
lines  of  the  flow,  thus  by  measuring  particle 
characteristics,  we  are  actually  measuring 
the  medium*s  characteristics.  This  is  easily 
achieved  by  use  of  particles  which  are  on  the 
order  of  a  few  microns  in  diameter. 

(2)  particle  generation   results   in  a  consistent 

particle  size,  or  at  least  a  majority  of  the 
particles  are  within  a  specified  range  of 
sizes.  For  example,  DOP  (di-2-ethylhexyl 
phthalate)  particle  generators  produce 
particles  with  a  mean  diameter  of  0.75 
microns  and  burning  turbine  oil  produces 
particles  on  the  order  of  4-6  microns  (from 
electron  microscope) . 


A  laser  is  incorporated  because  of  its  coherency  and 
power.  The  aforementioned  grid  is  in  reality  the  result  of 
two  coherent  teams  interfering  with  themselves  (Fig.  1) . 
The  optics  produce  the  beam  pair  by  conventional  splitting 
and  focusing.  Because  the  laser  beam  intensity  is  of  a 
Gaussian  distribution  the  grid  created  by  constructive  and 
destructive  interference  is  not  uniform  over  the  entire 
focal  volume.  The  fringes  created  on  the  outer  portions  of 
the  horizontal  axis  are  not  perfectly  formed  because  of 
slightly  different  intensities  in  that  area.  The  main 
portion  of  the  interference  grid  is  uniform  and,  most 
importantly,   the   fringe  spacing  (D  )  is  determined  only  by 

the  laser  wavelength  (lambda)  and  the  half  angle  of  the  beam 

pair  convergence  (theta/2) . 


D   =  lambda/ (2*sin (Theta/2)  ) 


The  number  of  fringes  (N  )  may  also  be  computed, 
provided  the  beam  focusing  convergence  angle  (delta  theta) 
is  known  ,by  the  relation: 


N  =8*tan (Theta/2) /(PI*Delta  Theta) 


The  number  of  fringes  is  important  only  to  ensure  that  the 
signal  processing  apparatus  will  receive  enough  information 
to  make  possible  further  data  reduction.  If  an  insufficient 
number  of  fringes  exists,  not  enough  signal  will  be 
available  for  velocity  determination. 

Solid  particles  with  a  mean  diameter  which  is  less  than 
half   the  fringe  spacing  are  introduced  into  the  free  stream 


FRINGE   PATTERN  GEOMETRY 


Figure  1 
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of  the  fluid  to  he  measured.  When  the  particle  is  small 
compared  to  the  fringe  spacing  the  scattered  light  exhibits 
essentially  100%  modulation,  while  particles  equal  to  the 
fringe  spacing  produce  negligible  modulation[ Ref .  2].  The 
particles  may  be  produced  by  atomization  of  water,  dop-air 
reaction,  burning  oil,  smoke  bombs,  or  may  be  bought 
commercially,  as  e.g.  polystyrene  spheres  and  some  solid 
chemicals  such  as  zinc  oxide.  Choice  of  the  particle  to  be 
used  is  normally  dictated  by  the  fringe  spacing  in  the  focal 
volume,  and  the  nature  of  the  medium  to  be  measured. 
Regardless  of  choice,  the  particle  must  be  small  enough  to 
follow  the  streamlines  as  closely  as  possible,  because  the 
laser  anemometer  actually  measures  the  particle  velocity  and 
unless  the  particle's  motion  is  a  true  representation  of  the 
fluid  velocity,  pertinent  data  may  not  be  obtained. 

Inside  the  measuring  (or  focal)  volume,  bands  of 
darkness  and  brightness  exist  due  to  the  interference  cf  the 
two  beams.  Particles  entering  the  focal  volume  will 
alternately  scatter  light  as  they  pass  through  regions  of 
constructive  interference  (  typical  focal  volume  dimensions 
are  300  microns  wide  by  4  millemeters  long  and  contain 
approximately  64  fringes) .  A  photomultiplier  tube  located 
on  the  optical  axis  and  focused  into  the  measuring  volume 
receives  this  scattered  light  and  converts  what  appears  to 
be  a  modulated  light  signal  into  a  series  of  voltage  pulses. 
Each  voltage  pulse  correlates  to  a  particle  crossing  one 
fringe.  Eecau-se  the  spacing  of  the  fringes  is  known,  the 
frequency  of  the  voltage  pulse  train  may  be  equated  to  the 
velocity  of  the  particle  by  the  following  relationship: 


f=V  /D   =  2*V   *sin(Theta/2) /lambda 
of      o 

where  D   is  the  fringe  spacing  and  V   is  the  velocity  of  the 
f  0 


particle.  Now  the  processor  is  left  with  a  fairly  simple 
task  of  counting  the  frequency  and  multiplying  it  by  a 
calibration  factor,  which  is,  in  essence,  a  fringe  spacing 
parameter  to  arrive  at  a  velocity. 

As  long  as  we  are  assured  that  only  one  particle  is  in 
the  measuring  volume  at  one  time,  the  apparatus  will  operate 
correctly,  but  this  is  an  unrealistic  assumption.  In  order 
to  achieve  such  a  situation,  the  seeding  density  of  the 
particles  would  be  exceptionally  sparse  and  the  resulting 
data  rate  would  be  so  slow  as  to  preclude  timely  data 
reduction.  To  overcome  this  obstacle,  a  validation  circuit 
is  incorporated  into  the  counter/processor.  When  a  particle 
crosses  a  fringe  spacing,  two  counters  are  triggered.  One 
measures  the  time  duration  for  five  fringe  crossings  and 
holds  this  information.  The  second  counter  continues  to 
measure  the  time  for  a  total  of  eight  fringe  crossings.  The 
time  of  the  first  counter  is  then  compared  with  that  of  the 
second  counter,  and  if  the  first  time  is  indeed  five-eights 
of  the  second  time  (within  specified  limits,  i.e.,  comp. 
accuracy)  a  validation  gate  is  opened  and  the  resultant 
frequency  is  then  stored  into  an  ensemble  cell.  The  ensemble 
cell  is  adjustable  to  hold  various  numbers  of  counts  and 
average  them  for  the  final  velocity  display.  The  final 
displayed  velocity  then  is  that  component  of  fluid  velocity 
which  is  perpendicular  to  the  optical  axis  and  contained  in 
the  same  plane  as  the  beam  pair  (Fig.  2)  . 

The  anemometer  may  be  operated  in  three  modes: 
reference  beam,  forward  scatter,  and  back  scatter.  The  last 
two  modes  are  most  frequently  employed  and  are  described 
here: 

Forward  scatter  mode:  in  this  mode,  the  photomultiplier 
tube  is  located  on  the  optical  axis  and  focused  in  a 
direction  directly  opposed  to  the  direction  of  the  emerging 
beam  pairs.    This   mode  offers   the   strongest   attainable 
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signal. 

Back  scatter  mode:  again  the  photomultiplier  tube  is 
mounted  on  the  optical  axis,  but  now  it  is  aligned  in  the 
direction  of  the  emerging  beam  pairs.  This  mode  affords  the 
ability  tc  measure  fluid  velocities  in  confined  spaces  or  in 
situations  where  the  dectector  cannot  be  mounted  in  forward 
scatter,  such  as  in  measuring  velocities  in  a  compressor. 
Additionally,  operation  in  this  mode  requires  a  higher 
powered  laser  due  to  the  reduced  intensity  of  back 
scattered  light  (Fig.  3)  . 
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Ill    ANALYTICAL  DERIVATION 

The  finite  element  method  is  a  numerical  analysis 
technique  fcr  obtaining  approximate  solutions  to  a  wide 
variety  of  enginering  problems.  It  offers  a  means  to  solve 
complex  continuum  problems  by  sub-dividing  them  into  simpler 
interrelated  problems.  The  nature  of  this  engineering 
problem  does  not  readily  lend  itself  to  a  closed  form 
solution.  Therefore  a  numerical  solution  is  dictated. 

The  method  of  finite  elements  is  well  employed  where 
there  are  irregular  geometries  or  unusual  specifications  of 
boundary  conditions.  The  method  gives  a  piecewice 
approximation  to  the  governing  equations.  Figure  4  shows 
the  finite  element  grid  employed-  in  this  flow  problem. 
Approximating  functions  (interpolation  functions)  are 
defined  in  terms  of  the  values  of  the  field  variables  at  the 
nodes.  For  the  finite  element  representation,  the  nodal 
values  of  the  field  variables  become  the  new  unknowns  and, 
once  found,  the  interpolation  functions  define  the  field 
variable  throughout  the  field  of  elements.  The  degree  of  the 
interpolation  polynomials  is  governed  by  the  number  of  nodes 
per  element  and  dictated  by  the  order  of  the  governing 
equations. 

The  basic  principle  of  the  finite  element  method  is  to 
divide  the  solution  domain  into  a  finite  number  of 
subdomains  (elements)  which  are  connected  only  at  node 
points.  Interpolation  polynomials  are  created  for  each 
element  and  a  collecton  of  interpolation  functions  over  the 
whole  solution  domain  provide  a  piecewise  approximation  to 
the  field  variable. 

Introduction  of  generalizations  facilitate  direct 
derivation  of  finite  element  equations  from  governing 
differential  equations.   The  method   of   weighted   residuals 


14 


•a 

•H 

■s 

n 

0> 

(D 

-P 

•d 

•H 


•H 


15 


(Gaierkii^s  Method)  is  a:  technique  for  obtaining  approximate 
solutions  to  non-linear  partial  differential  equations.  We 
first  assume  the  general  function  behavior  of  the  dependent 
field  variable  in  some  way  so  as  to  approximately  satisfy 
the  given  differential  equations  and  boundary  conditions. 
After  applying  this  approximation  to  the  original  equations, 
seme  error  will  result,  which  we  desire  to  minimize  or  drive 
to  zero  over  the  entire  solution  domain. 

The   basic   set  of  differential  equations  governing  the 

problem  from  a  velocity  and   pressure   formulation   approach 

are  (after  linearization  by  approximating  non-linear  viscous 

terms  u  ,  v  ) : 
n   n 


9u   3v    _ 
continuity   "^  +  g-  =  O 


momentum    u-n  3*  ■*■  vn  ~Z7.   s-pT  +  W  U. 


ay  "    p  9x 


u  ,  v  ,  and  E   are  approximate  solutions  to  the  flow  problem 
n   n       n 

th 
and  may  be  visualized  as  the  n    iteration  of  the   procedure 

of  solution.  These  equations  used  in  formulation  of  the 
problem  are  by  no  means  unique.  Stream  functions  and  stream 
function  with  vorticity  formulations  have  also  been  employed 
to  establish  finite  element  solution  procedures. 

Inasmuch   as   the  velocity  terms  are  of  second  order,  a 
second   order  interpolation   polynomial   must   be   used   to 
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calculate   these   unknowns.   The  pressure  terms  are  of  first 
order  and  hence  require  only  first  order  interpolation. 

The  basic  element  chosen  for  this  solution  was  the 
triangle  and,  in  a  localized  scheme,  is  visualized  as 
follows: 


Nodes   2,   4   and  6  lie  equidistant  between  their  respective 
corner  nodes. 

Nodes  1,  3,  and  5  are  referred  to  as  corner  nodes  and 
are  the  basis  for  the  linear  interpolation  polynomials 
derived  for  the  pressure  terms.  For  any  triangular  element, 
the   corner   nodes   are    described    by    their    Cartesian 

coordinates  x  ,  y  . 
i   i 

Consider  the  following  local  element: 


l — <y  >y  ) 


2    2 


<W 


The  coordinates  of  node  one  are  x  ,  v  ;   node   two   x  ,   y  ; 

11  2       2 

etc.   and   the   interior   point   (x,y)  may  be  described  as  a 
linear  combination  of  the  three  corner  nodes: 
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x=L  x  +L  X  +L  X 

11    2  2    3  3 


y=L  y  +L  y  +L  y 

11    2  2    3   3 


where    L    is   the   natural    coordinate   (interpolation 
i 

polynomial,  here  linear)  or  weighting  function  which  relates 

the   coordinates   of   the  corner  nodes  to  any  interior  node. 

L   has  a  value  of  one  at  the  point  x  ,y   and  is  zero  at   the 
l  i   i 

other  two  nodes  of  the  element. 


The  L  •  s  expressed  as  linear  functions  of  the  Cartesian 
i 


coordinates  then  are 


i 


Lt(x,y)=  ■%£   (ai+bix+ciy) 

L2(X'Y)=     lA     (a2+b22+C2y) 
L3(X'Y)=       IS     <VV+<V> 


where  2.A  = 


1  x  y 
11 

1  x   y 

2  2 

1  x  y 

3  3 


which  is  equal  to  the  area  of   the   triangular   element   and 
a  =x  y  -x  y  ,      b  =y  -y    ,  and  c  =y  -y  ,  the  determinants  of 

12332       123  112 

the  minors  corresponding  coefficients  a  ,  b  ,  or  c  ,  and  the 

ill 

other  coefficients  are  similarly  computed. 


For  the  velocity  terms,  second  order  interpolation 
polynomials  are  indicated,  while  the  pressure  terms  remain 
first  order.  Hence  the  introduction  of  more  nodes  (6)  will 
result  in  a  new  set  of  weighting  functions  which  may  also  be 
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expressed  in  terms  of  the  natural   coordinate   interpolation 
polynomials: 

N    =2L2-L  N    =4L    L 

111  212 

N    =2L*-L  N    =UL    L 

3  2         2  #23 

N    =2L2-L  N    =4L    L 

5  3         3  6  13 

P 

N      =L 

i       i 

e       4  e*  e       »     P 

u   =    £N   u  v  =    £N   v  P=2np       (xxxx) 

i     i   i  i     i   i  »  i 


Application   of   Galerkin's   Method   to   the   governing 
equations  results  in  the  following  intergral  expressions: 

here    Jft-         implies  "over  the  area  of  the   triangle". 


Ihen  we  may  make  the  substitution  of  the  field  variable 

in  terms  of  x  and  y  into  terms  of  L  . 

i 


e.g.:     given     u(x,y)    <=>   u(L  ) 

i 

then      §i±  -  tP— !  2!=J 

where     ^ti  r  _ki 

and  a  similar  procedure  for     dv  ~     ^***'   V 
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where    2!*  s  -£L 

e     e        e 
Substitution  of  equations  (xxxx)  into  u  ,   v   and   P   and 

application   of   the   above  chain  rule  allows  integration  by 
parts: 


v  j^fitt*  ^o-^e-vH^M 


K  = 

2 


0  9Nj 


-jNj^^y 


K  = 

3 


=  -Kly^y 


The  evaluated  integral  expressions  then  form  the  influence 
matrix  as  fellows: 


K 


-4— »■ 


K 


0  i  K 


K  '  K 


K  !  0 


'i   _J 


Specification  of  the  field  variables  or  right  hand  sides  of 
the  matrix  equation  then  provides  a  defined  problem: 

Additionally,  evaluation  of  the  integral  expressions  yields 
boundary  terms  which  must  be  computed  for  nodes  which  lie  on 
the  boundary  of  the  control  surface: 


.-(% 


B   = 


X  ds 


R   = 

2 


-J'i 


I  ds 


M 


where  X  =^7U'T1-Tixr      and   Y  =^W-Tl-nrr 
It  should  be  noted  that  along  vertical  boundaries   there   is 
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no  outward  pointing  normal  n   and  similarly  along  horizontal 

y 

boundaries  no  n   exists,  and  furthermore  these  terms  are  not 

x 

computed   where  u  and  v  are  specified.   Additionally,  it  was 

assumed  that  along  the  left-and  right-hand   side   and   upper 

boundaries  au/3*  and   /3y  were   identically   equal   to  zero 

because   these  boundaries   were   placed   at  five   or   more 

characteristic  lengths   (diameters)   from   the  cylinder  and 

considered  to  be  outside  the  realm  of  influence  of  the   flow 

about  the  cylinder. 

The  terms  comprising  K  ,  K   and  K   are  entered   element 

12         3 

by   element   into  the  main  program  in  steps  FLU02080  through 

QZ00170.  Then  each  element's  matrix  is  compiled  into  the 
field  matrix  (comprised  of  all  elemental  matricies)  and  is 
solved  iteratively  to  output  the  horizontal  velocity 
component  u,  the  vertical  component  v,  and  the  pressure  P  at 
every  node  in  the  field. 
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Hi     DESCRIPTION  OF  EXPE2IHENT 

A  right-angled,  three-sixteenths  inch  cylinder  made  of 
stainless  steel  was  mounted  in  an  AEROLAB  series  90  32x45 
inch  low  speed  wind  tunnel  through  the  floor  of  the  test 
section  (Fig.  5) .  The  portion  of  the  rod  protruding  down 
through  the  test  section  was  mounted  on  a 
two-degree-of -freedom,  micrometer  controlled  bed  which 
allowed  the  rod  to  be  maneuvered  in  the  vertical  plane 
parallel  to  the  flow  direction.  This  scheme  enabled  the 
laser  velocimeter  and  photomultiplier  tube  to  be  mounted, 
aligned  and  left  undisturbed  for  the  entire  period  during 
which  data  were  taken.  The  bed  was  positioned  to  read 
locations  in  the  control  volume  which  directly  coincided 
with  the  location  of  the  nodes  in  the  finite  element  mesh 
(Fig.  6)  . 

In  the  interest  of  precluding  as  many  opportunities  for 
error  as  possible,  all  u  velocity  components  were  first  read 
and  then  components  in  a  forty-five  degree  position  were 
taken  thus  necessitating  only  one  re-orientation  of  the  beam 
pairs. 

The  plexiglass  window  prefacing  the  optics  transducer 
of  the  test  section  had  one-quarter  inch  holes  drilled  into 
it  to  allow  the  beam  pairs  to  pass  unimpeded,  thereby 
minimizing  distortion  of  the  focal  volume  fringe  pattern. 
The  photomultiplier  tube  was  likewise  accommodated  through 
the  far  window. 

The  seeding  mechanism  employed  zinc-hexaclorcthane 
based  (non-toxic)  white  smoke  bombs.  The  nature  of  the 
closed  return  wind  tunnel  permitted  continuous  measurement 
without  interruption  for  re-seeding  of  the  flow.  By  means  of 
electron  microscopy,  the  average  size  of  the  smoke  particles 
was  determined  to  be  on  the  order  of  two  microns  (Fig.  7) . 
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Cylinder  mounted  in  wind  tunnel  and  LDV  setup. 


Figure  5 
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Traversing  mechanism 


Figure  6 


2^ 


Zinc-hexachlorathane  particles 
(X5000    lcm=  2,0/um) 


Figure  7 
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The  final  velocity  component  measured  at  each  node  was 
the  result  of  averaging  three  sucessive  readings  from  the 
velocimeter  processor.  Each  of  these  readings  represented 
the  average  of  4096  readings  which  had  been  validated  to 
within  1.5%  of  the  exact  five-eighths  validation  count 
criteria. 

The  flow  in  the  wind  tunnel  was  stabilized  at  five 
inches  per  second  to  maintain  a  Reynolds  number  of  about  40, 
which  is  in  the  realm  of  creeping  flow.  The  experimental 
setup  was  designed  to  measure  the  steady  flow  about  the 
cylinder,  hence  the  relatively  low  Reynolds  number,  and  no 
effort  was  made  to  obtain  unsteady  flow  measurements.  Three 
separate  runs  were  conducted,  two  over  the  forward  half  of 
the  control  volume  and  one  over  the  rear  half.  Because  the 
beams  of  the  laser  anemometer  were  approximately 
one-sixteenth  inch  in  diameter,  several  node  position 
measurements  were  excluded  due  to  inaccessibility.  The 
majority  of  these  were  in  close  proximity  to  the  cylinder. 
In  order  to  obtain  both  an  x-  and  a  y-velocity  component  of 
the  flow  at  every  node  which  was  accessible,  two 
measurements  were  made,  but  because  of  the  very  small 
vertical  component  present  in  most  of  the  flow,  direct 
vertical  beam  pair  measurement  was  unsuitable.  Instead, 
after  obtaining  the  horizontal  component  for  each  node  in  a 
run,  a  second  series  of  measurements  was  made  with  the  beam 
pairs  inclined  forty-five  degrees  from  the  vertical.  In  this 
fashion  the  y-component  was  determined  by  the  following 
relationship: 
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It  can  be  seen  here  that  the  unknown  resultant  (R)  has 
a  projection  onto  each  axis  representing  the  orientation  of 
the  beam  pairs,  thus  the  resultant  is  the  vector  which 
emanates  from  the  origin  and  terminates  at  the  intersection 
of  the  perpendiculars  constructed  from  the  u  and  v(at  45 
degrees)  vectors. 
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V.    DESCRIPTION  OF  COMPOTER  PROGRAM 

The  program  found  in  the  appendix  has  several  comment 
cards  inserted  for  clarification.  However,  a  broader 
presentation  is  offered  here.  The  program  may  be  subdivided 
into  four  main  sections: 

Section  1:  dimensioning  and  zeroing  of 
parameters  and  arrays.  (FL000020  through 
FLU00690) 

Section  2:   input  of  node  coordinates,  system 
topology  and  known  values  of   u,   v,   P,   and 
right   hand   sides,   and   output  of  the  same. 
(FLU00700  through  FLU02090) 

Section  3:  computation  of  interpolation 
polynomials  on  an  elemental  level  and 
arrangement  into  local  element  arrays. 
(FLU021C0  through  GO00320) 

Section  4:  loading  of  the  global  influence 
matrix  (TM)  and  solution  of  (TM)  (T)  =  (Q)  by 
Gauss-Jordan  partial  pivoting  and  back 
substitution,  with  an  iterative  routine 
incorporated  until  convergence  is  attained, 
and  associated  print  out.  (FLU04470  through 
FLU04880) 

The  computer  routine  systematically  tested  the  solution 
of  field  variables  with  the  previous  solution  until  all 
terms  in  the  T  and  T1  arrays  agreed  within  0.001  (ESPILA) . 
If  the  convergence  test  failed,  the  process  was  repeated 
from  section  3. 

Initially,  the  program  was  run  on  an  IBM  360  in  FORTRAN 
CLG  mode,  then  an  ensuing  run  was  completed  on  the  FORTRAN 
HCLG  compiler  and  a  time  reduction  of  two-thirds  was 
realized.  The  program,  as  given  in  the  appendix,  required 
approximately  375K   bytes  of  core  to  compile  and  675K  bytes 
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to  execute  and  output.  Substantial  space  savings  may  be 
realized  ty  utilizing  banded  matrix  storage  methods  vice  the 
method  employed  here,  but  due  to  time  limitations  this 
avenue  was  not  fully  explored. 

The  program  attempts  to  solve  simultaneously  for  all 
the  field  variables.  However,  introduction  of  zeros  on  the 
main  diagonal  of  the  TM  influence  matrix  resulted  in  a 
singular  matrix.  Another  approach  was  explored  where  the 
equation  cf  continuity 


dx       By 


was  replaced  ty 


2V  r        [   ax  ay  a^ay  J 


and  the  resulting  terms  were  added  to  the  TM  matrix  and 
their  respective  right  hand  sides.  The  only  field  variables 
affected  by  this  form  of  the  continuity  equation  were  the 
pressure  variables,  and  depending  on  whether  the  local 
corner  node  involved  in  the  local  pressure  equation  was  on 
the  boundary  or  interior,  different  right  hand  sides  had  to 
be  evaluated.  This  is  the  purpose  of  steps  PLU04460  through 
GO00320.  The  simultaneous  solution  technique  for  all  field 
variables  applied  to  this  new  set  of  equations  did  not 
converge. 

A  third,  unimplemented   technique   may   be   considered. 
Given  the  original  matrix  equations: 

(1)  K  0   +   0   +K  P   =R 

1  2        1 

(2)  0    +K  V   +K  P   =R 

13        2 

(3)  K  U   +K  V   +   0   =0 

2  3 

Solution  for  u  and  v  will  result  in 

-l       -l 

(1)     0=-K    K  P  +K    R 

v  '  12       11 
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-1  -1 

(2)     V=-K    K  P  +K    R 

13  12 

and  multiplication  by  K   and  K   yields 

2  3 

•  i  -  i 

K  &=-K    K    K  P  +  K  K    R 

2  2   12  2   11 

-1  -1 

K  V=-K  K    K  P  +  K  K    R 

3  3   13  3   12 

Summation   of   these   equations  allows  substitution  for  the 
third  matrix  equation: 


Tk  K    K   +  K  K    K 

i-21      2       31      3  J 


-1  _1    i  -1  .1 

P  =  KK    R+KK    R 

2   1      1       3   1      2 


The  matrix  operating  on   P   may   be   calculated   en  an 

algebraic  level  and  entered,  in  terms  of  a  ,  b   and  c  ,  into 

i   i      i 

the  program  in  order  to  solve  for   the   pressures   at   every 

corner   node.   The  calculated  P's  then  would  be  entered  into 

the  field  variable  array  and  the  partial   pivoting   routine, 

which   involves  only  those  terms  where  the  field  variable  is 

not  specified,  would  operate   primarily   on   the   first   two 

matrix   equations.    In   previous   computer  trials  where  all 

corner  node  pressures  were  specified,  although  not  known   to 

be   correct,   calculations   were   properly   performed.   This 

method  would  allow   the   actual   pressures   present   in  the 

control   surface  to  be  entered  into  the  computer  routine  for 

further  calculation  of  the  velocity  components. 
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ILz  MALYSIS    OF    RESULTS 

No   computer  routine  output  was  obtained  for  comparison 

of  results,  due  to  the  inability  of  the  program  to  converge. 

As   such,   comparison   of   the  two  methods  cannot  be  offered 

here.   The  non-convergence  of  the  program   was   investigated 

from   several  aspects:   validity   of   the  non-linear  terms, 

correctness  cf  the  loading  of  the   TM   matrix  on   both   the 

elemental   and   global   levels,    correctness   of   input, 

performance  cf  the   partial   pivoting   subroutines   and  the 

iteration   scheme;   all   without  success.   In  test  cases  the 

program  initially  solved  for  the  unknowns   within   magnitude 

accuracy   for   a   majority   of   terms.    However,   sucessive 

iterations  resulted  in  oscillating   solutions,   even   after 

twenty   iterations.    Convergence   was  expected  after  ten  to 

fifteen  iterations  without  under-  or  overrelaxation   factors 

incorporated   into   the   iteration   scheme.    After   the  new 

formulaticn  of  the  continuity  equation  was  incorporated  into 

the   computer   routine,  the  program  execution  again  resulted 

in  a  singular  matrix  when  only  the  left  and  right   hand   and 

upper    control  surface  boundaries   had   their   respective 

pressures  specified.   Underflow  occurred  in  the  execution  of 

the    subroutine    SOLVE    at    step   SOL00260.    Subsequent 

th         th 
investigation  revealed  that  the  178    and  179   rows  of   the 

reduced  matrix  equations  were  the  cause  of  the  underflow. 

The  data  which  were  taken  for  runs  number  two  and  three 
are  displayed  in  the  appendix  and  are  plotted  on  the  finite 
element  grid  in  Figure  8.  Without  the  analytical  results  to 
compare  the  measured  data  only  a  qualitative  statement  may 
be  made  pertaining  to  the  LDV  output.  It  is  readily  seen 
that  the  trend  of  flow  directions  around  the  cylinder 
appears  to  be  correct.  By  and  large,  the  flow  in  close 
proximity  tc  the  cylinder  was  accelerated  over  the  tcp  and 
decelerated  near  the  stagnation  points.   It  is   pointed   out 


31 


here  that  the  plot  is  a  composite  of  two  different  trials, 
and  as  such,  different  free  stream  velocities  were  present 
for  each  run.  In  Figure  8  the  vector  representing  the  flow 
at  a  node  was  plotted  and  the  magnitude  of  the  vector  in 
inches  per  second  was  displayed  above  the  vector  or  by  the 
node  from  which  the  velocity  vector  emanated. 
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A.      EREOB  ANALYSIS 

ae-positioning  of  the  cylinder  to  its  base  location  was 
accomplished  with  relatively  low  position  error.  Prior  to 
each  run  the  cylinder  was  aligned  using  the  beam  pairs  as 
reference.  In  the  horizontal  mode,  the  beam  pairs  would 
strike  the  leading  and  trailing  portions  of  the  cylinder  and 
reflect  to  positions  horizontally  opposed  on  the  far  side  of 
the  wind  tunnel  when  correctly  zeroed  vertically.  The  same 
procedure  was  repeated  for  the  vertical  beam  pair  except 
position  was  checked  by  the  upper  beam  reflecting  to  a 
vertical  position  on  the  far  wall  and  the  lower  beam  (at 
reduced  intensity)  striking  the  cylinder  support  on 
centerline.  This  procedure  was  repeated  several  times  and 
the  maximum  deviation  between  any  aligned  positions  was  0.2 
millimeters  horizontally  or  vertically.  Thus,  it  was 
determined  that  the  maximum  positioning  error  inherent  in 
the  system  would  be  the  ratio  of  the  smallest  coordinate  of 
any  node  measured  and  the  aforementioned  0.2  millimeters  or 

0.2/3.75=0.5% 

Freestream  velocity  measurements  were  likewise  tested 
for  fluctuations.  Typical  measurement  variations  ranged 
from  0.129  to  0.142  (run  number  three;  freestream 
measurements)  in  the  horizontal  mode  and  0.0947  to  0.100  in 
the  forty-  five  degree  mode  (velocities  in  meters/seconds) . 
These  figures  represent  a  total  of  10%  variation  which  may 
be  considered  as  ±  5.%  from  the  mean. 
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VII .  CONCLUSION 

The  experiment,  as  described,  required  approximately 
two  to  four  hours  for  each  run.  aside  from  the  initial 
alignment  of  the  components  and  the  cylinder,  this  was 
perhaps  the  greatest  drawback  of  the  system.  Next,  the 
exclusion  of  regions  in  close  proximity  to  the  cylinder  from 
measurement  was  deemed  to  be  a  major  factor  in  the  small 
amount  of  data  actually  taken.  The  latter  may  be  improved 
with  a  different  experimental  arrangement;  particularly,  the 
plexiglass  section  directly  prefacing  the  optics  transducer 
should  be  completely  removed,  thus  allowing  all  angles  of 
beam  pair  inclination,  and  in  turn  availing  more  locations 
to  measurement.  Another  improvement  which  could  be 
implemented  would  be  unitization  of  the  laser,  transducer 
and  photomultiplier  tube  arrangement.  This  would  decrease 
the  LDV ■ s  sensitivity  to  vibration  and  provide  a  more 
accurate  optical  axis  alignment.  The  aforementioned 
improvements  would  still  not  decrease  the  time  required  for 
data  collection,  but  the  introduction  of  a  two  component 
system  [Ref.4]  would  allow  simultaneous  measurement  using 
orthogonal  beam  pairs.  Furthermore,  this  type  of  system 
would  not  concern  itself  with  the  errors  introduced  in  beam 
pair  re-alignments  performed  with  the  single  beam  pair   LDV. 

The  reproducibility  of  the  experiment  was  judged  to  be 
good  with  the  exception  of  exact  duplication  of  the  wind 
tunnel  velocity  between  successive  trials.  This,  of  course, 
is  primarily  a  function  of  the  controllability  of  the  tunnel 
motor  and  propellor  assembly.  It  is  estimated  that  with 
adequate  facilities  for  incorporation  of  some  of  the 
aforementioned  improvements,  more  accurate  data  may  be 
obtained.  The  LDV  displays  great  potential  as  a  means  of 
measuring  a  wide  variety  of  dynamic  fluid  regimes  accurately 
and  without  introduction  of  disturbance  to  the  flow   itself. 
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APPENDIX  A 

LDV  DATA 

RUN  #2 

Node 
64 

U(cm/sec) 
12.43 

m45°( 
8.43 

cm/sec)     R( in/sec) 
4.90 

Z.0 

^73 

65 

12.33 

8.49 

4.85 

-1.5 

66 

13.10 

8.33 

5.18 

-5.8 

67 

13.17 

8.25 

5.22 

-6.5 

68 

13.87 

8.72 

5.49 

-6.3 

69 

13.87 

8.56 

5.50 

-7.2 

73 

14.03 

8.58 

5.57 

•7.7 

75 

13.40 

8.83 

5.29 

-3.9 

76 

13.43 

8.62 

5.31 

-5.3 

77 

13.00 

8.62 

5.13 

-3.6 

78 

13.47 

8.87 

5.32 

-4.5 

79 

12.87 

8.58 

5.08 

-3.3 

80 

13.63 

8.91 

5.38 

-4.3 

81 

13.20 

8.02 

5.25 

-8.0 

96 

10.50 

6.94 

4.14 

-3.7 

97 

12.83 

8.64 

5.06 

-2.7 

98 

13.47 

9.54 

5.30 

0.1 

99 

13.27 

8.92 

5.23 

-2.8 

100 

13.43 

9.25 

5.29 

-1.5 

101 

12.33 

8.83 

4.85 

0.7 

102 

9.77 

6.69 

3.85 

-1.8 

105 

13.03 

9.32 

5.12 

0.8 

106 

11.93 

8.68 

4.70 

1.7 
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RUN  #3 


Node 

U(cm/sec) 

V@4  5°  (cm/sec) 

R( in/sec) 

o 

1 

13.73 

10.77 

5.44 

o.23 

2 

13.70 

10.43 

5.43 

6.25 

3 

13.77 

10.33 

5.^3 

3.49 

k> 

13.23 

10.03 

5.22 

4.11 

6 

13.77 

10.93 

5.46 

7.00 

7 

13.80 

10.40 

5.44 

3.77 

8 

13.87 

10.17 

5.46 

2.11 

9 

13.27 

10.08 

5.24 

4.27 

13 

13.50 

10.47 

5.34 

5.54 

14 

13.57 

10.40 

5.36 

4.80 

15 

13.53 

10.20 

5.34 

3.76 

16 

12.30 

9.93 

4.89 

8.05 

19 

13.70 

10.43 

5.41 

4.38 

20 

14.60 

10.87 

5.76 

3.02 

21 

13.83 

10.57 

5.46 

4.63 

22 

13.60 

10.06 

5.36 

2.65 

23 

13.57 

10.23 

5.35 

3.79 

24 

12.93 

10.40 

5.14 

7.84 

25 

12.50 

9.92 

4.96 

7.00 

26 

11.90 

9.64 

4.74 

8.27 

29 

13.13 

10.16 

5.19 

5.40 

30 

12.57 

9.80 

4.98 

5.86 

31 

11.30 

9.76 

4.56 

12.50 

35 

14.17 

10.47 

5.58 

2.60 

36 

14.03 

10.43 

5.53 

4.36 

37 

13.23 

10.07 

5.22 

4.36 

38 

12.97 

10.09 

5.13 

5.72 

39 

11.83 

10.03 

4.75 

11.24 

5^ 

14.03 

10.87 

5.55 

5.58 

55 

14.17 

11.03 

5.61 

5.76 

56 

15.20 

10.57 

5.98 

-0.90 

57 

15.17 

10.47 

5.97 

-1.36 

58 

16.27 

10.00 

6.46 

-7.46 

59 

14.37 

10.20 

5.66 

0.20 
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